Singlet versus triplet particle-hole condensates in quantum oscillations in cuprates 
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Quantum oscillations in a tilted magnetic field offer the possibility of distinguishing singlet versus triplet 
order parameters in the particle hole channel provided the measurements reflect a putative "normal" state of a 
density wave obtained by applying a high magnetic field at low temperatures. A theoretical analysis is given 
that compares spin density wave, a singlet d-density wave, and a triplet d-density wave. While the existence 
of a spin zero in the oscillation amplitude is a necessary consequence of a singlet order parameter, a triplet 
order parameter may or may not exhibit a spin zero, making it a quantitative issue that depends on the actual 
extremal orbits on the Fermi surface. Nonetheless, a theoretical analysis can shed light on the striking recent 
measurements in YBa2Cua06+«. 



I. INTRODUCTION 

A remarkable experiment in 2007 1 involving quantum os- 
cillations of the Hall coefficient in YBa2Cu306.5 has raised 
an important question regarding the ground state, as the su- 
perconducting dome is crushed by a high magnetic field, high 
enough to destroy superconductivity to a large extent but not 
high enough on the scale of electronic energies. 2 3 Subse- 
quently a large number of such experiments have provided 
overall consistency, but perhaps not full agreement in all re- 
spects. Another striking observation has been a similar mea- 
surement in electron doped Nd2_ x Ce x Cu04 where the mea- 
surements are carried out between 30 — 65 Tesla, far above the 
upper critical field H C 2 less than 10 TeslaP- If any universal- 
ity is to hold, the mechanism of quantum oscillation must be 
the same in electron and hole doped cuprates, 5 even though 
the measurements in hole doped cuprates are perhaps below 
H C 2, but unquestionably carried out in the resistive state. An 
emerging view is that the oscillations result from Fermi pock- 
ets formed by a suitable density wave state, a condensate in the 
particle-hole channel. We explore this possibility to shed light 
on two conflicting experiments in a tilted magnetic field, one 
in which no spin zeros were found 7 and the other in which 
they were found P 

Condensates in particle-hole channel are fundamentally dif- 
ferent from condensates in the particle-particle channel — 
superconductors. Since there is no exchange requirement be- 
tween a particle and a hole, the symmetry of the orbital wave 
function does not constrain the symmetry of the spin wave 
function. A given orbital symmetry can come in both singlet 
and triplet varieties. These order parameters are very ni cely 
classified on the basis of the angular momentum channelPSl 
Here, we concentrate on the s and d-orbital channels leaving 
out the possible p-channel, which probably is unlikely for the 
cuprates. Thus, the singlet (sDDW) and the triplet (tDDW) d- 
density wave order parameters are, considering only the two- 
dimensional and the two-fold commensurate case: 

(ct(k+Q) C/3 (k)) = i$ Q f(k)S a0 , (1) 
( C t(k+Q) C/J (k)) = i$ Q /(k)a • <r aP . (2) 

Here c Q (k) is the fermion destruction operator of spin index 



a and the d-wave form factor is 

/(k) = (cos k x a — cos k y a) /2. (3) 

The vector Q = (w/a, ir/a) and the magnitude $q is real. 
The singlet order parameter transforms as identity in the spin 
space, while the triplet transforms as h ■ er Q( g, n being the 
direction of the spin quantization axis and cr the standard Pauli 
matrix. We shall assume that spin-orbit coupling is negligibly 
small. 

The singlet and the triplet s-wave density waves are the 
conventional charge (CDW) and spin density waves (SDW), 
defined by setting /(k) = 1 and removing the factor of i: 

(4(k + Q) C/3 (k)) = $ Q< 5 Q/3 , (4) 
(4(k + Q) C/3 (k)) = $qn-^. (5) 

Note that we have denoted the magnitudes of all the order 
parameters by $q for notational simplicity. 

The sDDW modulates neither charge nor spin but repre- 
sents staggered circulating charge currents, as shown in Fig.[T| 
In this case, the order parameter breaks lattice translational 




FIG. 1 : (Color online) Circulating charge current for a single square 
plaquette corresponding to singlet DDW order parameter. Note that 
there is no circulating spin current. 



symmetry, time reversal, parity, and a rotation by ir/2, but the 
product of any two symmetry operations is preserved. The 
order parameter is "hidden" because most external probes do 
not couple to currents or variations of kinetic energy. The or- 
der parameter corresponding to S = 0, as in a sDDW, does 
not respond to magnetic field directly, as in a paramagnet, but 
the quasiparticles energies are split by the Zeeman effect. 
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The tDDW is further hidden because it is invariant under 
time reversal unlike SDW. Nonetheless both break spin rota- 
tional symmetry and lead to Goldstone modes. The tDDW, as 
well as the SDW, undergo a spin-flop transition in an arbitrary 
small magnetic field for zero spin-orbit interaction. If there 
is anisotropy in the spin space from spin-orbit coupling, there 
will be a non-zero threshold field beyond which the spin-flop 
transition will take place. In cuprates spin-orbit coupling is 
small enough that for high magnetic fields relevant for quan- 
tum oscillations, it is almost certain that the spins will be per- 
pendicular to the applied magnetic field H, as shown in Fig. [2] 
There is a particularly nice way to bring out the similarities of 
tDDW and SDW, if we define the macroscopic order parame- 
ters by 

y = i^/(k)^ -^<4(k + Q)c^(k)}, (6) 

k a,P 

m = ^5>^(4(k + Q) C/3 (k)). (7) 

k a,0 

Note the trace operation in the above definitions. The tDDW 
and the spin flop in the presence of a magnetic field is shown 
in Fig. 
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FIG. 2: (Color online) Spin-flop of triplet DDW in a magnetic field 
H pictured for a single square CuO plaquette with O atoms ommited. 
There is circulating spin current but no charge current. Here (a) cor- 
responds to the situation when H = and (b) to H ^= 0. 



n. SPIN ZEROS 

In this section we shall calculate the spin zeros in quan- 
tum oscillations following the Lifhsitz-Kosevich (LK) formal- 



ism|2] which we do not duplicate here. Only the relevant as- 
pects pertaining to the reconstructed Fermi surfaces due to 
sDDW, tDDW and SDW will be described. In addition, we re- 
strict ourselves to two dimensions without c-axis warping and 
bilayer splitting in YBa2Cti306+,5.^This simplicity reveals 
the essential elements without unnecessary clutter. The exten- 
sion to include the neglected effects are entirely straightfor- 
ward. However, we do need to specify the various parameters 
needed for the illustrative calculations, and these are summa- 
rized in the Appendix|A] In this appendix we also provide the 
frequencies in the absence of the Zeeman field. 



A. Singlet 

The spectra of excitations for sDDW can be obtained from 
the following Hamiltonian 



UeRBZ 



(8) 



where the reduced Brillouin zone (RBZ) is bounded by k y ± 



±7r/a and 



e t)k iW k 

-iW k e Tik+Q 

e |)k iW k 

-iW w e |>k+Q 



(9) 



f,k = £k + § VbH and e lM = e k - \[i B H, where \i B 



Here e 

is the Bohr magneton. The four component spinor is given by 



L k+Q,t' C k,V 



,t 

'k+CUV 



Of course the up and the 



down spin sectors are decoupled, and the eigenvalues are 



\ s — 
a„ l- — 



- ek- 
2 



QJ 



± 



EkH 

4 



QJ 



<) 



(10) 



which clearly shows Zeeman splitting as for a free spin, 
what follows, to simplify the notation, we shall define: 

(\ = 2~[ e k — Ek+Cj] 



C2 



= \ [ek + ek+cj] 



In 

(11) 
(12) 



and drop the wave vector arguments when there is not any 
possibility of confusion. For notational simplicity, we denote 
a = +1 =t and a = -1 =|. 

The Zeeman term depends on the total magnetic field. It is 
included exactly within the Hamiltonian, because it is crucial 
in determining the sensitive interference of the amplitudes of 
quantum oscillations related to spin. The effect of the mag- 
netic field on the orbital part is treated differently. First, its 
effect on the DDW order parameter can be neglected, as long 
as the system is deep inside the DDW phase; the fields rele- 
vant to quantum oscillations are energetically weak perturba- 
tions on the DDW gap. 13 Clearly, close to a quantum phase 
transition where the DDW gap collapses, this will no longer 
be true — a situation that is not relevant to the experiments ad- 
dressed here. Second, in calculating quantum oscillations us- 
ing LK formalism only the extremal Fermi surface areas and 
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its various derivatives are needed, which can be calculated in 
the absence of the orbital part of the magnetic field. The oscil- 
lations arising form the extremal in-plane orbits are of course 
determined by only the normal component H z ; for simplicity 
we are ignoring warping and bilayer splitting. At the wave 
vector (7r/2a, 7r/2a), and at symmetry related points in the 
Brillouin zone, there are nodes. If the chemical potential, /i, 
is placed at these nodes, the orbital diamagnetic susceptibility 
diverges as H — > in the absence of any interlayer coupling 
and at zero temperature, T = 0. 10 However, the values of \x 
considered here are too far from the nodal points for this to be 
relevant. 

From the LK formalism it is immediately obvious that the 
spin interference factor, R s , within the LK formalism is 



COS Tip 



m* g 1 



m 2 cos 9 



(13) 



where cos 9 is the angle between the magnetic field with re- 
spect to the normal and p stands for the p-th harmonic. Within 
our mean field theory g = 2. The singlet DDW order pa- 
rameter must exhibit spin zeros corresponding to elementary 
fermionic excitations of charge e, spin 1/2, g — 2, and the 
cyclotron mass m* calculated from Eq. [10] Nonetheless, the 
system is not a paramagnet because the spectra in Eq. [^cor- 
respond to a broken symmetry state far from a conventional 
paramagnet but with total spin 5 = 0. Within Fermi liquid 
theory, with the reconstructed Fermi surface, Eq. 13 should 
hold to all orders in perturbation theory. However, residual 
Fermi liquid interactions between the quasiparticles can cer- 
tainly lead to renormalization of m* and g. However, quite 
generally, within the Fermi liquid formalism, the residual in- 
teractions can only increase m* and g, assuming small spin 
orbit coupling, which appears to be the case; see Ref. [8] and 
reference therein. Thus the case for the existence of spin zeros 
will be stronger if the Fermi liquid corrections are taken into 
account. 



do not have definite spin unlike the singlet case. For arbitrary 
k the mixing depends on the k-space orbit, and this gives rise 
to a dynamically generated spin-orbit effect, as can be seen 
from the quasiparticle spectra, if we expand in powers of the 
magnetic field. The differences between the energies of the 
two close energy levels, either the holes or the electrons, are 
given by 



A, 



A,. 



gfi B H- 



(16) 



V^! + W 2 + {g^ B H/2f 
g eS (k)» B H, (17) 



that is, by an effective g-factor. A typical plot of g e s is shown 
in Fig. [3] Another way to picture the distinction between the 




FIG. 3: (Color online) A map of g e ff(k) for tDDW. Note that it is 
mostly 2 for most of the hole pocket and mostly zero for the electron 
pocket. This effective g-factor will split the hole and the electron 
pockets into two close pockets depending on the magnetic field. The 
actual splitting due to Zeeman coupling is not visible on the scale of 
the plot. 



B. Triplet 

In the triplet case the spin orientation is chosen to be per- 
pendicular to the direction of the applied field because of the 
spin flop. Thus, with no loss of generality, we can choose 

H z = H cos 9, H x = H y = and n x = 1, n y = n z = 0. 
Then the matrix A is 



A = 



The eigenvalues are 











iW 





e t.k+Q 


-iW 








iW 







-iW 








e i,k+Q 



(14) 



e 2 ± ^{ tl + ^cj^ B HY + W 2 



(15) 



The four eigenvectors are various linear superpositions involv- 
ing the coherence factors of the original fermion operators and 



singlet and the triplet cases is the schematic energy level dia- 
gram shown in Fig. [4] 



1. Effective g-factor 

Let us focus on the electron pockets; the results for the hole 
pockets follow identically. Let 



£± = £2 



£l ± |M-B-ff 



w 



1/2 



(18) 



To calculate the interference factor due to the two close elec- 
tron pockets within the LK formalism, we must compute the 
sum 



3=± 



(19) 



4 
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FIG. 4: The splitting of the levels in a magnetic field. Note that in the 
triplet case we cannot identify the levels with a definite spin quantum 
number. The hole and the electron pockets are labelled "hole" and 
"elec". 



corresponding to the Landau level index n m ± corresponding 
to the extremal Fermi surface orbits. Since the splitting is 
small for experimentally relevant H, it is a good approxima- 
tion to assume that n m+ (/i) « n m _ (/.i) w n m (/i). Simi- 
larly, 



dn r 



de 



+ 



( dn r - 



\ 9e_ 



where n m = 
e — Ae. Thus, 



(n m+ + n m -)/2, e + 



27ripn m (/i) 



<9e 



(20) 



e + Ae, and e_ 



x cos 27rp 




(21) 



We can set (e — /i) w and the identification of e and Ae 



follows trivially from Eq. 18 Then the interference factor is 

to* 5eff(k) 1 



R s = cos irp 



cos 



(22) 



As the <7 e ff is k-dependent we average over the extremal orbits 
shown in Fig. [3] 



2. Filling of the Landau levels 

There is another method to obtain the interference factor 
that serves as a consistency check, especially because we 
have made approximations. In this approach we focus on 



the filling of the Landau levels. We can define the average 
n = I (n+ + n_) and the difference An = | (n+ — n_). 
Where the plus sign makes reference to the pocket that is 
bigger and the minus sign to the pocket that is smaller. It 
does not matter if the difference is due to the Zeeman spin 
splitting, as for the singlet case, or the splitting in the spec- 
tra in the triplet case. If we extend this notation to the ex- 
tremal areas to, n m = ^ (n m + + n m _ ) and the difference 



h (n r , 



-) Again taking advantage of small 



splitting, we can show that in the LK formula 



(23) 



3=± 



Now the dependences on the g e ff, cyclotron masses, and the tilt 
angle are contained in An m . If we compare with the previous 
results we obtained for the cosine factor we get 



cos np 



9s,t m 



1 



2 TOp COS ( 



= cos (2irpAn m ) . (24) 



valid for both singlet and triplet. If An m = the ar- 
gument vanishes and the cosine is unity; there will be no 
spin zeros at all. We can deduce what is the effective ce- 
faclor for the triplet if we know the result for the singlet: 
(An m ) t / (An m ) s = g t m* t j g s m* s . Assuming that the effec- 
tive cyclotron masses are approximately the same (see Table[l]) 
for sDDW and tDDW, we get 



(An m ) t 
(An m ) s 



(25) 



TABLE I: Calculated cyclotron masses in units of the free electron 
mass for sDDW, tDDW, and SDW order parameters, e-p and h-p 
correspond respectively to electron and hole pockets, and the two 
rows to the two distinct Fermi surfaces. The errors are less than one 
part in 1000. The first row corresponds to the inner pocket and the 
second to the outer pocket. 



sDDW 
e-p 


sDDW 
h-p 


tDDW 
e-p 


tDDW 
h-p 


SDW 
e-p 


SDW 
h-p 


1.55 
1.55 


0.95 
0.97 


1.56 
1.54 


0.96 
0.96 


1.27 
1.28 


1.12 
1.13 



It is easy to see that that n m is proportional to the extremal 
Fermi surface areas. For our purposes we can set g s = 2 to 
compute g t from the integration of the Fermi surface areas. 
The results are shown below in Table |Hl The absolute value of 
R s , the factor that is responsible for the spin zeros is shown 
for the electron pocket, as a function of the tilt of the magnetic 
field with respect to the normal in Fig. [5] In the singlet case 
we have chosen the free electron value of g s = 2. In the triplet 
case the calculated effective g t = 0.32 is used for the electron 
pocket. We can see that for the singlet a spin zero is found 
close to 50°, while the first spin zero for the triplet occurs 
above 60°. Similar results are shown for the hole pocket in 
Fig. fusing g t = 1.65. 
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TABLE II: The effective g-factors. The computations were carried 
out for 40 Tesla field for illustrative purposes. The Method 1 is based 
on <? e ff(k) in Eq. 22 and the Method 2 is based on Eq. '. 



25 



A reason- 
able estimate of the errors may be taken to be the difference of the 
results between the two methods. 



3eff 


tDDW 
e-p 


tDDW 
h-p 


SDW 
e-p 


SDW 
h-p 


Method 1 
Method 2 


0.32 
0.31 


1.65 
1.65 


0.45 
0.40 


1.42 
1.42 
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FIG. 5: (Color online) \R S \ as a function of the angle of the tilted 
magnetic field for the electron pocket and for p — 1. The solid line 
corresponds to sDDW and the dashed the tDDW. 



C. Commensurate SDW 

A commensurate SDW in our scheme is an orbital s-wave 
order parameter for which the calculation proceeds identically 
except that the /(k) = 1. The results are shown in Fig. [7] 
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FIG. 6: (Color online)|i? s | as a function of the angle of the tilted 
magnetic field for the hole pocket and p — l.The solid line corre- 
sponds to sDDW and the dashed the tDDW. 



FIG. 7: (Color online) \R a \ as a function of the angle of the tilted 
magnetic field for SDW and p = 1. The solid line in this figure 
represents the hole pocket and the dashed line the electron pocket. 



III. CONCLUSION 



The rigorous statement we can make is that a singlet order 



parameter will definitely lead to Zeeman splitting (see Eq. 10 1 



and therefore spin zeros (recall the discussion of the g-factor 
above). A triplet order parameter, be it SDW or triplet-DDW, 



does not exhibit Zeeman splitting (see Eq. 15 i. It does lead to 
some dependence of the amplitude of oscillations on the tilt 
in a non-universal manner, and may or may not lead to spin 
zeros within a reasonable angular range, such that the system 
is not driven into the superconducting state, which may lead 
to the loss of the amplitude of the oscillations. In other words, 
the experiments in Ref. [8] are consistent with a singlet order 
parameter. 

It is remarkable that once the conventional band parameters 
are adopted, and the oscillation frequency is adjusted close 
to the experimentally observed one at 526 T for the electron 
pocket (for the other frequency, effects such as bilayer split- 
ting and c-axis warping, must be taken into account), the os- 
cillation amplitude as a function of the tilt angle for the sin- 
glet order parameter is very close to the measured amplitude 
in Ref. 8 As one can see from the plots, the triplet order pa- 
rameter behaves very differently. 

The spin zeros determine only the product gm* /m e , which 
calculated in the case of sDDW for the electron pocket is 3.1, 
as opposed to 3.2 in Ref.|8] For the hole pocket, in contrast, it 
is calculated to be 1.92. Unfortunately there does not appear 
to be any evidence of the hole pockets in experiments, about 
which we have commented elsewhere!^ 

A useful insight into tDDW can be obtained from the co- 
herence factors corresponding to the four eigenvectors of the 



matrix in Eq. 14 These are 



71, k 

72, k 

73, k 

74, k 



Vi,kCk,t + "i,kCk+cu, 

«2,kCk+Q,t + u 2,kCk,i, 
«3,kCk+Q,t + U3,kCk,i, 
^4,kCk,t + M 4,kCk+Q,i, 



6 



where the coherence factors are 



and 




Ae kja = ei + ^a/igH. (29) 

The coherence factors reflect the fact that the quasiparticles 
do not have a definite spin. For arbitrary k the mixing is not 
of equal amplitude. 

Another amusing observation is that as far as the product 
gm* I m e is concerned both tDDW and SDW yield essentially 
identical answers even though the individual values for the ef- 
fective g-factors and the m* are different. For the electron 
pocket this product is 0.5 for tDDW and 0.51 for SDW. Sim- 
ilarly, for the hole pocket it is 1.59 for tDDW and 1.60 for 
SDW. 

One might wonder if there are other possibilities of a singlet 
order parameter that could be a candidate broken symmetry 
state. In principle one cannot rule out CDW order, although 
such an order parameter would not have gone unnoticed in 
many other direct experiments. Thus, the existence of spin 
zeros may be consistent with sDDW, whose direct observation 
by its very nature may be considerably hidden. 
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Note added: After our work was completed we learned of 
two interesting papers that appear to be complementary. In 
one, the problem has been addressed from the perspective of 
stripe physics (arXiv: 1007.1047vl) and in the other antifer- 
romagnets are studied (arXiv:1006.0167vl). 



Appendix A: Parameters 

We use a common band stmctureP^l 

£k = — 2t(cos k x a + cos k y a) + At' cos k x a cos k y a 

— 2i"(cos 2k x a + cos 2k y a), 
where t = 0.154 eV, t' = 0.32<, t" = 0.5*'. For sDDW and 
tDDW we choos d 9 ' 10 ' 17 ' 



Wk = —z~ (cos k x a — cos k y a), (A2) 



with Wq = 0.9t. The chemical potential ji = — 0.775t then 
results in a hole doping of approximately 10.7%. For SDW the 
choices were Wk = Wo — 0.675t and /i = — 0.93t, yielding a 
hole doping of approximately 10.6%. The frequencies corre- 
sponding to the Fermi surface areas in the absence of Zeeman 
splitting are (a) sDDW: 513 Tesla for the electron pocket and 
1005 Tesla for the hole pocket; (b) tDDW: 513 Tesla for the 
electron pocket and 1004 Tesla for the hole pocket; (c) SDW: 
539 Tesla for electron pocket and 1012 for the hole pocket. 
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